Particle distribution and nuclear stopping in Au-Au collisions at v /J/v r /v=200 GeV 



L.L.Zhu and C.B. Yang 
Institute of Particle Physics, Hua-Zhong Normal University, Wuhan 430079, P.R. China 

The transverse momentum distribution of produced charged particles is investigated for gold-gold 
collisions at ^/sjvjv = 200 GeV. A simple parameterization is suggested for the particle distribution 
based on the nuclear stopping effect. The model can fit very well both the transverse momentum 
distributions at different pseudo-rapidities and the pseudo-rapidity distributions at different cen- 
tralities. The ratio of rapidity distributions for peripheral and central collisions is calculated and 
compared with the data. 
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I. INTRODUCTION 

Ultra-relativistic heavy ion collision is the only way to 
study nuclear matters at extremely high temperature and 
density in the controlled experimental conditions. Such a 
study can provide information about the internal struc- 
ture of hadrons and the interactions among partons from 
which those hadrons are formed. An important quantity 
for characterizing the properties of the produced parti- 
cles is the transverse momentum distribution in different 
rapidity regions. The transverse momentum distribution 
at mid-rapidity region can tell us the matter density cre- 
ated in heavy ion collisions and that at large rapidity tells 
us the baryon stopping strength in the collisions. Experi- 
mentally, the charged hadrons can be comparatively eas- 
ily identified and investigated. The multiplicity distribu- 
tion of charged particle is often used as a global measure 
of the dynamics of the ultra-relativistic heavy ion colli- 
sions. The total multiplicity as a function of rapidity or 
pseudo-rapidity can be used to measure the nuclear stop- 
ping effect in the collisions. Because the stopping effect 
is directly related to how effectively the kinetic energy 
of the colliding particles can be converted into thermal 
one and how many secondaries can be produced in the 
process, it has been a focus on both nucleon-nuclear and 
nucleus-nucleus collisions for quite a long time and been 
studied in various approaches before 0, 12, IE El IE IE • 
Some of the investigations are done at the nucleon level, 
some others are based on the string model. Since it is 
not feasible to perform a first-principle parameter-free 
calculation of the nuclear stopping effect, experimental 
guidance plays a crucially important role in such studies. 

This paper discusses the transverse momentum dis- 
tributions of emitted particle at the highest energy of 
the Relativistic High Ion Collider (RHIC) between gold- 
gold (Au-Au) at ^/snn — 200 GeV for different central- 
ities and rapidities. We will start from a very simple 
parameterization for the particle distribution stimulated 
from momentum degradation effect and fit the experi- 
mental data of the transverse momentum distributions 
for charged particles for different centralities and rapidi- 
ties. Then the rapidity distributions can be obtained 
from integration over transverse momentum for different 
centralities and are compared to the experimental data. 



The organization of the paper is as follows. In sec- 
tion II, we will propose a simple parameterization for the 
invariant particle distributions. Then in section III, we 
fit the experimental data to our model parameterization. 
The last section is for a brief conclusion. 



II. PARAMETERIZATION FOR THE 
INVARIANT PARTICLE DISTRIBUTION 

More than 30 years ago, the transverse momentum dis- 
tribution of particles produced at polar angle — ir/2 in 
high energy proton-proton collisions was suggested || as 

E^=Aexp(-e/2b) , (I) 

from the analogy of Landau's hydrodynamical 
model, with £ = \n(mT + Pt)/™ the transverse ra- 
pidity and b a parameter to be determined by fitting the 
experimental data . Here m is the mass of the produced 
particle, tut = \J m 2 + p\ is the transverse mass. Re- 
cently this parameterization was also used in describing 
the data for nucleus- nucleus collisions at RHIC energy 
in the mid-rapidity region. At small px the above 
transverse momentum distribution can be approximated 
by an exponential and the collective flow effect can 
be extracted if the temperature of phase transition is 
known. So, The value of b can tell us some information 
about the collective flow. In recent years, some authors 
have used the blast wave model to parameterize the 
transverse momentum distribution, see |10| for example. 
That model is valid only at low px region, since the 
blast wave is a collective dynamical phenomenon for 
the produced hot medium and the contribution from 
hard partons to the particle production cannot be 
contained in the blast wave model. For the distributions 
of particles with nonzero longitudinal momentum px , no 
first-principle theory, such as quantum chromodynamics, 
can be used to derive a formula for the distribution. 
Therefore, a phenomenological ansatz is needed. In 
this paper, as usually done in other studies, we simply 
assume a factorized form for the distribution as 

rf 3 N 

E—=Aexp(-e/2b)H( PL ). (2) 
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In a former study of p— A collisions it is shown that 
the Feynman xp (or longitudinal momentum p^) distri- 
bution can be approximated by an exponential form. The 
exponential decrease of the longitudinal momentum dis- 
tribution is a reflection of the momentum degradation 
effect from nucleon-nucleon interactions while incident 
nucleon is penetrating the target. For nucleus- nucleus 
collisions, one can expect that the distribution may be 
not too far from that for p — A collisions. Thus we as- 
sume 



H(p L ) <x exp(-cpf) , 



(3) 



with c, d two parameters characterizing the longitudinal 
momentum degradation effect in nucleus-nucleus colli- 
sions. For p— A collisions d = 1 If the physics for the 
longitudinal momentum degradation is the same in p — A 
and A — A collisions, d ~ 1 can be expected. A small 
deviation from d = 1 in A — A collisions may be due to 
the multiple NN collisions, each of which has a different 
center of mass energy. Here, parameter c is a measure 
of the strength of longitudinal momentum degradation 
effect. 

By changing variable from rapidity y to pseudo- 
rapidity i], using p l = pr sinh(?y), one can get the particle 
distribution as follows 



dN 



exp( 



1'Kp r pdp r pdTj 
ln 2 (22±£T)) 

v m 11 

2b 



Apr cosh(?7) 
n x + Pt sinh 2 (^) 

c(pt sinh(?7)) d ) . 



(4) 



We assume that the four parameters A, b, c and d de- 
pend only on centrality of the collisions at given collision 
energy, but not on the pseudo-rapidity. 
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FIG. 1: The charged particle distribution dN/2irpTdpTdri as 
a function of transverse momentum pr for Au+Au collisions 
at pseudo-rapidity rj — (top panel) and rj — 2.2 (bottom 
panel) for four centralities. The data points are taken from 
|12jl . The solid curves are the fitted results with Eq. For 
clarity, the spectra have multiplied by the indicated factors. 



III. CENTRALITY DEPENDENCE OF 
PARTICLE DISTRIBUTIONS 

To determine the parameters for different centralities, 
we fit the experimental data for produced particles with 
different centralities and rapidities. For this purpose 
we choose the data from BRAHMS Collaboration for 
Au+Au collisions at ^/smn — 200 GeV 12]. It needs 
to be mentioned that Eq. 0} is supposed to be valid 
for a specific species of hadrons with given mass m. In 
the chosen data, only the spectra of charged particles 
(h + + h~)/2 are given, where all charged particles with 
different masses are included. If the contribution to the 
spectra can be parameterized by Eq. (0} for each species 
of charged hadrons, the distribution for charged particles 
should be described by the sum of a few terms, each of 
which is given by Eq. Q . Then there would be too many 
parameters. Instead, in this paper, we assume that Eq. 
Q is still valid for all charged particles, as long as m in 
Eq. is taken as an effective mass. In principle, the 
effective mass should depend on the relative contents of 



all charged particles. It is known that the proton to pion 
ratio is larger than 1 at px ~ 3GeV/c [l3| and that the 
kaon to pion ratio is about 0.2 in a wide range of pr- So 
the contribution to the distribution of charged particles 
from protons and kaons will be quite large, especially for 
central collisions in the mid-rapidity region. Thus the rel- 
ative contents may be different for collisions at different 
centralities. Even with the same centrality, the relative 
contents change with px and r/. So generally speaking, 
the value of m is a function of centrality, pt and r\. For 
simplicity, we assume in this paper that the value of m 
in Eq. 0} is a constant for all centralities, px and rj. 
For this reason, we choose, for different centralities and 
rapidities, the mass m in Eq. 0} to be 0.51 GeV, about 
the average mass of proton, kaon and pion. With the 
charged particle px distributions at rj = and 2.2 at the 
same centrality cut, the fitting process can be completed 
in two steps. We first fit the experimental distributions 
of transverse momentum at mid-rapidity (ry = 0) to de- 
termine the parameters A and b, then, with fixed A and 
b, get parameters c and d from fitting transverse momen- 
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centrality 


A 


b 


c 


d 


~ 10% 


239 


0.351 


0.0878 


0.94 


10 - 20% 


163 


0.350 


0.0826 


0.94 


20 ~ 40% 


93.3 


0.344 


0.0818 


0.94 


40 ~ 60% 


39.8 


0.340 


0.0788 


0.94 



TABLE I: Values of parameters A, b, c and d for distribution 
of dN /2nprdpTdri for four centrality cuts. The fitted data 
are from 113. 



turn distributions at T) = 2.2 at the same centrality. The 
experimental data chosen have a pt range from 0.45 to 
4.65 GeV/c. We found from fitting that the value of d 
is 0.94, very close to 1 as expected, and almost the same 
for different centralities. So we fix d — 0.94 for the dis- 
tributions at other centrality cuts. The fitted results are 
shown in Fig.l together with the experimental data from 
BRAHMS Collaboration for pseudo-rapidity 77=0 and 2.2 
at four centralities. It is obvious that the data can be fit- 
ted very well to the simple ansatz shown in Eq. @. The 
parameters obtained are tabulated in TABLE I. 

In TABLE I the value of b decreases slowly from cen- 
tral to peripheral collisions. This means that Rcp = 
{dN / NcPrdpT)c I (dN / NcPrdpr) p increases with pt at 
t] = 0. For pions, Rcp decreases with pt for pr not too 
small, while Rcp for protons increases up to px about 
4 GeV/c p. 4]. The behavior of Rcp for unidentified 
hadrons is somewhat in between those for pions and pro- 
tons. The parameter c is an indicator of nuclear stopping 
in Au-Au collisions. The value of c dictates the differ- 
ence between the spectra at different pseudo-rapidities at 
given centrality. The larger the value of c, the bigger the 
difference, thus the stronger the nuclear stopping effect. 
One sees that the nuclear stopping effect gets weaker for 
peripheral collisions. This is in agreement with the naive 
expectation that the incident nucleons need to traverse 
more nucleons in central collisions and lose more energy. 
Though the values of c are quite small, the little differ- 
ence can give the difference of the px distributions at 
different pseudo-rapidity rj. 

After obtaining the values of parameters A, b, c and d 
in Eq. Q for the particle distribution, we can get the 
pseudo-rapidity distribution for different centralities 



dN 
dr) 



p T dp 7 



dN 



PTdprdrj 



(5) 



where the factor 2 is included for the total number of 
all charged particles. BRAHMS Collaboration have pre- 
sented the pseudo-rapidity distributions for 6 centralities 
[l5|. Unfortunately, the centrality cuts in are differ- 
ent from those in [l2| . So the obtained parameters shown 
in TABLE I cannot be used for a comparison of our the- 
oretical expectation and the data on pseudo-rapidity dis- 
tributions. Furthermore, we would try to calculate the 
ratio between rapidity distributions for peripheral and 
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FIG. 2: Fitted invariant charged particle yields as a function 
of pseudo-rapidity r\ for six different centralities in Au+Au 
collisions at y^ijvjv = 200 GeV. The data points are taken 
from HHil . Solid curves are from our fitting with Eq. 



centrality 


A 


b 


c 


d 


0-6% 


513.3 


0.540 


0.21 


0.94 


6 ~ 15% 


387.9 


0.540 


0.19 


0.94 


15 ~ 25% 


268.6 


0.541 


0.18 


0.94 


25 ~ 35% 


188.6 


0.541 


0.175 


0.94 


35 ~ 45% 


119.4 


0.542 


0.16 


0.94 


45 ~ 55% 


69.4 


0.542 


0.154 


0.94 



TABLE II: Values of parameters A, b, c and d for distribution 
of dN/drj for six centrality ranges. The fitted data arc taken 
from lla|. 



central collisions and compare it with experimental data 
which is given in 0] by PHOBOS Collaboration. So we 
now choose to fit the PHOBOS data on rapidity distri- 
bution | lGj . One should keep in mind that the pseudo- 
rapidity distribution of charged particles includes con- 
tributions from charged particles with all possible trans- 
verse momentum for any specific pseudo-rapidity. So soft 
particles with very low px contribute most to r\ distri- 
butions. Since the multiplicity of pions is always much 
larger than those for protons and kaons when the trans- 
verse momentum is integrated, pions always dominate 
the pseudo-rapidity distribution. For this reason, the 
mass m in Eq. (JHJ) should be very close to that for pions. 
So we take m = = 0.139 GeV in the calculation of 
the pseudo-rapidity. Because the value of m is different 
in fitting the px and rj distributions, other parameters 
will also be different. However, we will fix d = 0.94 in 
fitting the rj distributions. The data on pseudo-rapidity 
distribution for the charged particle are over the pseudo- 
rapidity range —5.4 < 77 < 5.4 for six different centrality 
cuts. The fitted curves are shown in Fig. 2 together with 
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the data points. The agreement is excellent. The fitted 
parameters are tabulated in TABLE II. The value of b is 
almost constant for six centralities, but c decreases faster 
from central to peripheral collisions than shown in TA- 
BLE I. The small increase of parameter b from central 
to peripheral collisions is responsible for the decrease of 
the suppression effect for pions in the collisions. Such a 
trend is opposite to that in TABLE I, because in Fig. 1 
protons contribute more in central collisions than in pe- 
ripheral ones and the production of proton is, instead of 
being suppressed, enhanced in central collisions. When 
we focus on the pseudo-rapidity distributions, the effect 
of c is to narrow the distributions. In our model such 
effect is mainly from nuclear stopping. The fact that the 
value of c in TABLE II (where protons play less impor- 
tant role) is much larger than that in TABLE I suggests 
that the nucleus stopping effect can be illustrated more 
easily from pion spectra than from proton spectra. 



is measured in PHOBOS experiment at different colliding 
energies . With the fitted parameters for the pseudo- 
rapidity distributions, the ratio Rpc can easily be calcu- 
lated. We use the parameters for centrality cut 35 — 45% 
in the calculation of Rpc, with N p given in |16|. and 
the obtained results are plotted in Fig. 3 in comparison 
with the data points with the centrality cut 35 — 40% 
for the peripheral collisions. Up to rf — our calculated 
results agree with the data quite well. For rj > 0, only 
data points at ^/sWn = 19.6 GeV are available. In col- 
lisions at such a low energy, the contamination from the 
spectator nucleons to the produced particles in the very 
forward region is considerably serious, especially for pe- 
ripheral collisions. This contamination will make Rpc 
larger. That is a possible origin of the discrepancy of our 
calculated results from the data at large rj '. To verify 
this statement, more deliberate experimental investiga- 
tions are required. 
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FIG. 3: The ratio Rpc of dN c h/dr]' per participant pair be- 
tween peripheral (35-40%) and central (0-6%) collisions. The 
data are taken from [l6|. The solid curve is for the calculated 
results from parameterization Eq. @ and Eq. JSJ. 



IV. CONCLUSION 

In summary, we have investigated the transverse mo- 
mentum and pseudo-rapidity dependence of the distri- 
bution of charged particle multiply for several centrality 
ranges, respectively. We assumed a factorized parame- 
terization for the invariant momentum distribution, with 
the transverse momentum distribution at mid-rapidity 
being given by a Gaussian form of the transverse rapidity 
and the longitudinal momentum distribution by a expo- 
nential form characterizing the nuclear stopping effect. 
For the produced charged particles, with suitably chosen 
effective mass, the simple phenomenological factorized 
form can fit the experimental pt and rj distributions very 
well. Then the enhancement of the rapidity distribution 
for peripheral collisions relative to central ones can be 
explained naturally for very forward particle production. 



One can go further to the very forward region with 
pseudo-rapidity ij ~ j/bcam = 5.36 for ^/snn — 200 
GeV Au+Au collisions to investigate the centrality de- 
pendence of rapidity distributions. To compare the dis- 
tributions at different colliding energies, a shifted vari- 
able is introduced rj' = r; — j/beam? and an observable 

R ( 1 dN\ I ( 1 dN\ 

\N P dr) J 35 _ 40% / \N P d v J _ 6% 
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